Absence of isentropic expansion in various inflation models 
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Dynamics of the inflaton is studied when it interacts with boson and fermion fields and in minimal 
supersymmetric models. This encompasses multifield inflation models, such as hybrid inflation, and 
typical reheating models. For much of the parameter space conducive to inflation, the inflaton is 
found to dissipate adequate radiation to have observational effects on density perturbations and 
in cases to significantly affect inflaton evolution. Thus, for many inflation models, believed to 
yield exclusively isentropic inflation, the parameter space now divides into regimes of isentropic and 
nonisentropic inflation. 
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I. INTRODUCTION 



The standard picture of inflation 0, 0, Hj] has two distinct regimes. In the first phase, the universe undergoes 
inflationary expansion at constant entropy, so meaning with no radiation production. Thus if inflation started at 
finite temperature, such as in new inflation scenarios 0, the universe would supercool during the inflation phase, and 
£Ni ■ if inflation started near zero temperature, such as in chaotic inflation scenarios |3j, the universe would remain there 
during the inflation phase. Then in the second phase of inflation, radiation is reintroduced or introduced in a period 
called reheating. In this letter we show that for a large number of commonly studied inflation models, this picture is 
flawed. We show that in many models significant radiation production occurs during the inflation period, which has 
noticeable observational effects on energy density fluctuations as well as in cases completely altering the evolution 

■ history of the inflaton. 

Oh! We examine the model of a scalar inflaton field in interaction with other scalar \ and fermion fields through the 

■ basic interaction terms 



where hereafter cf> is the inflaton field, \ and V'x are additional fields to which the inflaton couples, and tpd are light 
fermions into which the scalar v-particles can decay rru, > 2m^, d . Aside from the last term, this is the typical 
Lagrangian used in studies of reheating after inflation 0, 0, or more to the point, if the inflaton is to eventually 
release its energy into the universe, such interactions are to be expected. Moreover in multifield inflation models, such 
as hybrid inflation fj| , the bosonic interaction plays a basic role in defining the model. 

We can qualitatively divide the study of the interactions into two regimes, weak coupling g,g' < 10~ 4 and 
moderate to large perturbative coupling 10~ 4 < g,g' < 1. In the weak coupling regime the radiative corrections are 
too small to affect the flatness of the the inflaton effective potential. In the moderate to large perturbative coupling 
regime, the radiative corrections from individual terms are significant, so some mechanism, like supersymmetry, must 
be invoked to cancel these effects and maintain an ultraflat inflaton effective potential. Later in the paper we will 
show that in minimal SUSY extensions of the typical reheating model, decay channels for the x 01 ^x particles are 
present and the ipd field above is simply a representative example. Since in the moderate to strong perturbative 
regime, reheating models will require SUSY for controlling radiative corrections, Eq. with inclusion of the ipd field 
thus is a toy model representative of the typical reheating model. 

Before proceeding, the regime in which the radiation energy density during inflation is below the Hubble scale 
p r < H 4 will be referred to as isentropic, supercooled,or cold inflation [l], 0> Ll El EJ and the other regime p r > H A 
will be referred to as nonisentropic or warm inflation Q • As referred to Eq. , this letter will focus on the moderate 
to large perturbative coupling regime. 

In the conventional approach 0, 0, 0, 0, 0; f° r the effective evolution equation of the inflaton background 
component, (p = {(p), the assumption has been that aside from radiative corrections that modify the ^-effective 
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potential, V e s(ip), this equation is the same as its classical counterpart. However, we have shown in earlier works 
HUGS, t na t i n addition to radiative corrections, quantum effects also arise in the (^-effective equation of motion 
(EOM) in terms of temporally nonlocal terms, which generically lead to dissipative effects. Moreover, although SUSY 
cancels large quantum effects in the local limit, for the dynamical problem the nonlocal quantum effects can not be 
cancelled by SUSY. Our previous considerations of inflaton dynamics were limited to nonexpanding spacetime. Here 
we have extended the calculation to the expanding case (for related earlier works see 0>E3)- We also provide in this 
work possible physically viable particle physics models, including those based on SUSY, that realize the dissipative 
dynamics being considered. In this letter we will not delve too much into the (sometimes) evolved technical details 
of the calculations and the quantum field theory dynamics, for which we refer the interested reader to our previous 
papers 01 an d the most recent one 13]. Here we will focus on the application and relevance of the calculations to 
inflation and cosmology in general. For this, in Sec. [H]we start by giving the results for the inflaton effective equation 
of motion that emerges from a response theory derivation that generalizes to expanding spacetime our previous 
calculations [T^j. Then in Sec. IIHI we present physically viable particle physics models that can realize the typical 
interactions and decaying mechanism leading to our results, discussing in particular models based on supersymmetry 
where large quantum corrections, that otherwise could spoils the usual flatness requirements for the inflaton potential, 
can be kept under control. A minimal SUSY model is explicitly presented and its properties discussed. Finally the 
conclusions are given in Sec. IIVI 



II. THE EQUATION OF MOTION AND DISSIPATION IN DE SITTER SPACETIME 

Our calculation here is a generalization to expanding spacetime of our previous results flU ] and whose details are 
also reported in [l^. The starting point is the full equation of motion for the inflaton, taken as a homogeneous 
classical background, ip = <p(t), which for the couplings shown in Eq. gives 

Cp(t) + 3H(t)<p(t) + £R(*Mf) + + <7V(*XX 2 > + 9'(Tp x i>x) = , (2) 

where H{t) = a/a is the Hubble parameter, R(£) the curvature scalar, and £ the gravitational coupling. V(p) is the 
tree level potential for the inflaton. In the calculations that follow, we will choose this to be a quartic potential A0 4 /4. 
To satisfy density perturbation constraints, it is well known the self coupling is tiny, A ~ 10 -13 , which is why in 
writing Eq. J5J we have neglected the insignificant quantum fluctuation corrections coming from the self-interaction 
of the inflaton field. Thus for our demonstration of dissipation and radiation production, in the ip effective EOM 
it is enough to consider the inflaton field as a classical background field interacting with a scalar x an d spinor ip x 
quantum field. As shown previously [HI, the two quantum correction terms in Eq. J5J lead to terms contributing 
both in the linear regime with respect to the inflaton amplitude from the Yukawa interaction due to the fermionic 
quantum corrections and in the nonlinear regime from the quadratic coupling of \ to the inflaton. 

So far we have not said anything about the role of the ipd spinors in Eq. Q and their contribution in Eq. J2J. 
Their main effect is to dress the propagator for the \ scalar and so they enter nonperturbatively in Eq. J2J through 
the loops of the x field. In addition, if we also allow open decay channels of \ into ipd, 4>di with decay rate r x _^ j; , 
the main effect of these spinors is to provide an effective damping in the dressed x Green's functions. This will in 
turn reflect back in the ip effective EOM Eq. as an effective dissipation for the inflaton field arising from the 
g 2 p(t)(x 2 } term. If the inflaton is also allowed to decay into the i\) x fermions, this can manifest itself in Eq. J2J as 
a dissipation term due to the inflaton direct decay and will be important for instance in the linear regime for the 
inflaton amplitude, or more specifically, at the time of reheating when the inflaton may oscillate around the minimum 
of its effective potential. These latter effects that can come from the fermionic quantum correction in Eq. J2J are not 
the effects we arc interested in here; these effects have been well described in several previous studies of reheating 
after inflation 0,0- Here we are concerned instead with the dissipation identified in Refs. 0,^1, associated with 
the nonlinear regime for the inflaton field, which can manifest during inflation; these effects mainly arise from the 
nonlocal quantum corrections due to the x scalars in Eq. Pjl. As shown in our previous references, we can then 
express Eq. @ in terms of a nonlocal effective EOM, relevant in the nonlinear regime, in the form 



(pit) + 3H(t)<p(t)+£R(t)<p(t) 



dV cS (p(t)) 
dtp(t) 



+ Ac/pit) / dt'<p(t')<p(t')K(t,t') = 0, 

J t 



(3) 
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where ^3| 

»«W = [qV(t)/« 2 (T) + m\ + 2(6? - lJH 1 ] 1 ' 2 , (5) 

m x — gip ^> m^ d , and r X)t (r) ~ /i 2 rn^/[87rw Xit (r)] is the decay width of scalars x of comoving momentum q into 
fermions ip^. The scale factor has been chosen a(t) = exp(Ht), since we restrict the analysis to the inflationary, 
quasi-deSitter regime. H = \J SvrV'off / (3m^j) is the Hubble parameter, with mp\ the Planck mass, and V e g((p(t)) is 
the renormalized effective potential for the inflaton, after the local quantum corrections to V(ip) have been taken into 
account. 

The kernel K(t,t'), Eq. J3J, is obtained by a response theory approximation, which takes advantage of the slow 
dynamics of the inflaton field so that the dynamics happen effectively in the adiabatic regime. This is equivalent to 
the closed time path formalism at leading nontrivial order 8, 9] , which treats the effect of the field x on the evolution 
of <p(t). The mode functions for the x _ field that enter in the Green's functions appearing in the response theory 
expansion of the (x 2 ) term in Eq. (J2J) (see e.g. Ref. [H3) are computed from a WKB approximation, which treats 
expansion effects and the time variation of the background inflaton field. These mode functions are then used in 
the loop calculations that determine the effect of the x quantum corrections to the background inflaton EOM. This 
WKB approximation is generally valid in the adiabatic regime doj x /dr <C w 2 . In our calculation this adiabatic regime 
emerges since for the parameter values we will be studying the ip motion is overdamped and since the x mass is heavy 
m x ^> H. In fact, these two conditions, the slow ip dynamics and the heavy x mass are enough to assure the validity 
of a WKB approximation for the x modes. In addition, the condition m x 3> H implies that curvature effects in the x 
field quantum corrections to the background inflaton field are subleading, with dominant term being the Minkowski 
like corrections. This will be explicitly observed below when we numerically compare several results based on Eq. 
Note, in the limit H — > 0, a — ► constant, Eq. Q agrees with the corresponding kernel for nonexpanding spacetime 
m 1QJ. The physical origin of the nonlocal (dissipative) term in Eq. JSJ is as follows. The evolving background field 
ip changes the mass of the \ boson which results in the mixing of its positive and negative frequency modes. This in 
turn leads to coherent production of x particles which then decohere through decay into the lighter ?/v-fermions. 

When the motion of tp is slow, an adiabatic-Markovian approximation can be applied that converts Eq. (|3J 
to one that is completely local in time, albeit with time derivative terms (for Minkowski spacetime details of this 
approximation are in [lOj , while for the expanding spacetime case this is derived in 13] ) . The Markovian approximation 
amounts to substituting t' — ► t in the arguments of the (/^-fields in the nonlocal term in Eq. ©. The Markovian 
approximation then requires sclf-consistently that all macroscopic motion is slow on the scale of microscopic motion, 
thus (p/tp,H < T x . Moreover when H < m x , the kernel K(t,t') is well approximated by the nonexpanding limit 
H — > 0. The validity of all these approximations will be examined below. Combining these approximations, the 
effective EOM Eq. J3J becomes 



druJ x .t{T) 



exp 



-2/;„ drT x , t (r) 



(4) 



v5+[3tf + TM]^ + £R^+^M = 0, (6) 

dip 



where, by defining a x = h 2 /(Sir) 



Tfo) = ^ g4Q ^ 2 (7) 



647rm x ^/l + a 2 yyl + a 2 +l 

Fig. ^compares the various approximations for a representative case where g = h = 0.37 and the inflaton potential 
is that for chaotic inflation V e g(ip) = Ac/? 4 /4 with A = 10 -13 |3J. In Fig ^ evolution has been examined at the 
final stages of chaotic inflation where we start with <p(t = 0) = mpj. The solid line is the exact result based on 
numerically solving Eq. J2J). Plotted alongside this, although almost indiscernible, is the same solution expect using 
the nonexpanding spacetime kernel (dashed line), obtained by setting H — > 0, a — > constant in Eq. Q}, and the 
solution based on the adiabatic-Markovian approximation of Eq. © (dot-dashed line) for the same parameter set. 
As seen, the expanding and nonexpanding cases differ by very little and the adiabatic-Markovian approximation is in 



4 




FIG. 1: Evolution of ip(t) for A = 1CT 13 , g = h = 0.37, £ = 0, <p(0) = m P1 , £(0) = 0. 



good agreement with the exact solution. This confirms simplifying approximations claimed in |3, lift lll| but up to 
now had not been numerically verified. 

More interestingly, and the first major result of this letter, the dotted line in Fig. \T\is the solution that would be 
found by the conventional approach in which the nonlocal term in Eq. (JSJ is ignored. The conventional approach 
Qi IE 01111 El expects the inflaton to start oscillating, which is the precursor to entering various stages of pre/re- heating. 
However with account for dissipative effects, this never happens for our example in Fig. 1, since the inflaton remains 
overdamped till the end when it settles at its minima at <p = 0. Moreover, throughout inflation, and not just near the 
end, the inflaton dissipates energy, which yields a radiation component of magnitude 



Pr « Tp 2 /(4H) . (8) 

For the case in Fig. ^ the overdamped regime for the inflaton persists until g < 0.35, below which its evolution 
at the end of inflation has oscillatory features similar to conventional expectations. However the inflaton is still 
dissipating radiation at the level Eq. (JSJ) all throughout the inflation period. The radiation energy density produced 
through this dissipative mechanism is greater than the Hubble scale, p r > H A and T x > H, if for example g = h and 
g > 10~ 2 or, as another example, if h = 0.1 and g > 10 -3 . If this radiation thermalized, which is expected for T x > H, 
then T ~ p X J A > H and it is known that the primordial density perturbations, Sp/p, produced by the inflaton are 
altered from their zero temperature value |7 Jljj. Il5l Il6j- In particular in the standard relation Sp/p ~ HStp/tp, for 
T <H, Sip 2 ~ H 2 is the cold inflation result |fH W for T > H and T < 3H, Sip 2 ~ HT and for T, T/3 > H, 

Sip 2 ~ V HTT 16]. Applying these relations to obtain the spectral index n s for the A</> 4 /4 model, for cold inflation 
the known value is 1 — n s = 3/N e [TtI ]. whereas in contrast for strong dissipative warm inflation (T > 3H), for the 
same normalization wc find 1 — n s = 3/(2N e ), where N e is the number of e- folds of inflation. For example for N e = 60 
the discrepancy in n s between the warm and cold inflation cases is about 3%, thus possibly resolvable with current 
generation high-precision satellite experiments, WMAP and Planck. Moreover cold inflation for any e-folding N e > 1 
has the well known r/ and large ^-amplitude problems, since observationally consistent slow-roll solutions require 
?7i0 < H and p > mp\ respectively [17$. On the other hand, for warm inflation these problems are nonexistent, 
since we find for example under the observational constraints Sp/p = 10 -5 and N e = 60, m 2 ^ w 190-ff 2 > H 2 and 
(p ~ 0.087mpi < mp\. We thus arrive at our second major result. In multifield inflation models, or in fact any inflation 
model once reheating interactions are accounted for, there are parameter regimes feasible to inflation, that have never 
been properly treated since the nonlocal term in Eq. JJJl is neglected. Once this term is retained, it is seen that up 
to fairly weak coupling, adequate radiation is produced during inflation to alter density perturbations. Although this 
conclusion is based on situations where thermalization is assumed, one can infer the same qualitative outcome once 
p r > H 4 , irrespective of the statistical state. 

In regimes where Eq. (0 produces significant radiation during inflation, models such as in 0, 0, 0] have to be 
examined on a case by case basis to determine the effect on density perturbations. Dissipative effects themselves can 
produce a rich variety of spectra ranging between red and blue |lj, 0, |2(j . A dramatic example is in [2(j which 
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treats new inflation type symmetry breaking potentials. In cold inflation, such models yield relatively featureless red 
spectra, whereas with account for dissipative effects, the spectra can be altered between red and blue producing, as 
an example, the running blue to red spectra suggested by WMAP/2df/Sloan pH ]. 

III. SUSY MODEL IMPLEMENTATIONS AND PARTICLE DECAYING CHANNELS 

The above results rely on the ability of the x-particles to decay, which in our toy model Eq. is fulfilled by the 
i/jd-fermions. Based on Eq. (0J, this dissipative mechanism is not specific to what the x particles decay into, but 
simply requires a nonzero T x > H. Thus the typical possibilities for a scalar field are either decaying into fermions, as 
in our toy model, or into gauge bosons. Also, in models where x develops a nonzero vacuum expectation value, such 
as hybrid inflation the 4>~X interaction generates the term g{x)x ( t )2 1 which implies a direct x ~ * 00 decay channel. 
Independent of inflation issues, particle physics models are adept with interaction configurations similar to Eq. 
For example, the simplest implementation of the Higgs mechanism in the Standard Model has the background Higgs 
field coupled to W and Z bosons, thereby generating their masses, and these bosons then are coupled to light fermions 
through the well known charged and neutral current interactions. In the context of inflation, conventional studies 
typically only present the inflaton sector itself, since the tacit assumption in these studies has been that interactions 
are unimportant to the inflation phase. The lesson learned by this work is interactions can play a vital role within the 
inflation phase, and so care must be taken to understand how the inflaton sector is embedded within the full theory. 
The few studies that have attempted to do this Tgjcontain interaction structures similar to Eq. (Q. 

Whether in our model Eq. Q , reheating models [5| , or any other inflation model which has a cf> — x coupling in the 
moderate to large perturbative regime @, sizable corrections are induced to the (^-effective potential and so SUSY 
must be invoked. A minimal SUSY extension that incorporates the 4> — x coupling would be 

W= V\^> 3 + g<PX 2 + fX 3 + mX 2 , (9) 

where $ = <\> + t/jd + 9 2 F and X = x + Oi/j x + 2 F X are chiral superfields. In the above model, there would be no 
additional fermion to associate with ipd from our toy model Eq. However the x-ficld has a decay channel via a 
i/j x triangle- loop into two light inflaton bosons <fi. For this case, everything in Eqs. ©-I® is unaltered except the 
decay channel is different with now T x ~ (fg 2 ) 2 m x . Thus there are additional factors of coupling constants, but in 
the large perturbative regime, the evolution of ip should be similar to the solid line in Fig. ^ Similarly the ^> x -fermion 
can decay into and ip, mediated through a triangle loop comprised of heavy x and ipx particles. 

Other weaker channels also occur. For example, the x and ip x SUSY partners can have masses differing by Am > m^, 
without harming the one-loop cancellations in the ip effective potential. Thus in the model Eq.@, a x particle could 
decay into two fermions, like in our toy model Eq. with one being the fermion partner ^ x and the other an 
inflatino ip. This channel is phase space suppressed since m x — ra^ x ~ m^, which leads to T x ~ g 2 m^/rn x . Based on 
Eqs. © and JSJ, this channel leads to T > H for g > 0.3. Weaker channels such as these need not dominate radiation 
production but also can serve other particle physics functions such as baryogenesis 0, l23j or dark matter production. 

Up to now very little has been said about the role of the ^ x -fermion in our toy model. Here its only purpose is 
to supply a one loop contribution that can suppress the large contributions from the x sector (then four x fields are 
required in Eq. QJ), thus mimic SUSY For de Sitter spacetime, these results are modified due to ^-dependent 

mass corrections and the ko integration in this geometry. Both modifications add corrections from the Minkowski 
spacetime effective potential by terms ~ 0(<C l)g 2 H 2 ip 2 < m^ip 2 and so can be neglected. Additional analysis of 
both quantum and thermal corrections for SUSY models like Eq. has been studied independently by the authors 
of Ref . |24j . They have also concluded that quantum corrections can be kept under control and in addition they find 
the same holds for thermal corrections. 

As briefly discussed below Eq. ©, for warm inflation in the strong dissipative regime, an important result is 
the avoidance of the ry-problem. In cold inflation models, the ^-problem is that slow-roll inflation requires < H 
whereas particle physics effects and the coupling to the background metric often prohibit this from happening. In 
particular in SUSY models, since they imply supergravity, it means higher dimensional operators inherently enter 
in forms such as a n (f> n /m^ -4 or V<p 2 /m 2 , which at either large or small field amplitude can cause the 77-problem. 
However for warm inflation in the strong d issip ative regime, slow-roll and density perturbation requirements can be 
satisfied for > H and even 3> H 0, Il6l fl9L |20| , thus intrinsically there is no ^-problem in any model. 
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IV. CONCLUSIONS 



It is now opportune to give a few remarks on the applications and limitations of our calculation. Eqs. © and 
are valid for any decay channel for the %-boson, thus any T x , including the case of no decay channel, T x = 0. If 
T x > H, then the approximation Eq. © also applies. In the regime, T x < H, within the characteristic oscillation 
period of the free inflaton field ~ l/m^, K(t, t') oscillates considerably, since m x 3> m$, H and without attenuation, 
since T x < H. Numerically this is a more difficult regime for studying the net effect of the nonlocal term. Another 
point to note is that Eq. only has the temperature independent contribution to dissipation and so only provides a 
lower bound. For T x > H thermalization can in principle be justified and when the temperature is bigger than either 
the inflaton mass, T > m$, or the decay product mass, T > m^ d in our toy model, dissipative effects are enhanced 
pllfil2fij. In further work, all these points need to be investigated. Although this work focused on inflaton dissipative 
effects induced by a x-boson field, similar considerations can be applied for a ip x fermion. This implies that fcrmionic 
and bosonic pre/reheating may be equally affected by our results. In any event, for moderate to large perturbative 
coupling, SUSY forces both types of interactions to be present side-by-side. 

Inflation and reheating phases have been in the past treated as almost detached problems, as so has the quantum- 
to-classical and the initial condition problems and so on. On the other hand, the warm inflation picture shows, here 
and elsewhere [?|,E3>II3j that inflation can be treated as a single consolidated scenario. As one outcome, various 
problems, namely 77, graceful exit, quantum-to-classical transition, large inflaton amplitude, and aspects of initial 
conditions |2(1 |2jj , can be nonexistent in warm inflation. 
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